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FREE VIBRATION OF UNIFORM AND NON-UNIFORM EULER BEAMS USING 
THE DIFFERENTIAL TRANSFORMATION METHOD 


MAHMOUD A. A., ABDELGHANY S. M., EWIS K. M. 


Abstract. In this paper, the differential transformation method (DTM) is applied for free vibration analysis of 
beams with uniform and non-uniform cross sections. Natural frequencies and corresponding normalized 
mode shapes are calculated for different cases of cross section and boundary conditions. MATLAB code is 
designed to solve the differential equation of the beam using the differential transformation method. 
Comparison of the present results with the previous solutions proves the accuracy and versatility of the 
presented paper. 


1. Introduction 


In order to calculate fundamental natural frequencies and the corresponding mode 
shapes, different variational techniques such as Rayleigh_Ritz and Galerkin methods had 
been applied in the past. Besides these techniques, another numerical methods were also 
successfully applied to beam vibration analysis such as finite element method. 


The DTM is a technique that uses Taylor series for the solution of differential 
equations in the form of a polynomial. The Taylor series method is computationally 
tedious for high order equations. The differential transform method leads to an iterative 
procedure for obtaining an analytic series solutions of functional equations. Since the 
beginning of 1986, Zhou and Pukhov [1] have developed a so-called differential 
transformation method (DTM) for electrical circuits problems. In recent years researchers 
had applied the method to various linear and nonlinear problems. Ayaz [2] applied The 
differential transform method to the differential algebraic equations , Moustafa El-Shahed 
[3] applied DTM to the non-linear oscillatory systems , Reza Attarnejad and Ahmad Shahba 
[4] applied DTM to the free vibration analysis of rotating non-prismatic beams, 
Rajasekaran S. [5] used the differential transformation (DT) to determine the natural 
frequency of beams and columns, Reza Attarnejad et al. [6] presented an application of 
differential transform in free vibration analysis of Timoshenko beams resting on two- 
parameter elastic foundation. Differential transform method is applied to the sine-Gordon 
equation by Jafar Biazar and Fatemeh Mohammadi [7], which arise in differential geometry, 
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propagation of magnetic flux and stability of fluid motions. J. Biazar and M. Eslami [8] used 
DTM for Nonlinear Parabolic-hyperbolic Partial differential equations. DTM was applied to 
linear and nonlinear system of ordinary differential equations by Farshid Mirzaee [9]. 
Keivan Torabi et al. [10] applied DTM for longitudinal vibration analysis of beams with 


non-uniform cross section. The non-linear vibration analysis of beams is applied by Qiang 
Guo and Hongzhi Zhong [11] using a spline-based differential quadrate method. 


In this paper, the vibration problems of uniform and non-uniform Euler-Bernoulli 
beams have been solved analytically using DTM for various end conditions. 
2. Basic Idea of Differential Transformation Method 


Following ref [2] we can obtain the idea of the DTM, The differential transformation of 
function u (x)is defined as follows; 





(1) 


dx* 


_ 1( d‘u(x) 
a ! 


In Eq. (1), u(x) is the original function and U(k) is the transformed function. Differential 
inverse transform of U/(k) is defined as follows; 


u(x) = yx'UW) (2) 


k=0 


In fact, from (1) and (2), we obtain 





u(x) “52 ( eo) (3) 


Eq. (3) implies that the concept of differential transformation is derived from the Taylor 
series expansion. From the definitions (1) and (2), it is easy to obtain the following 
mathematical operations; [2, 3, 4, 6, 9] 


1.If f(x) = g(x) +h(x), then F(k)=G(k)+H(k). 
2. If f(x)=cg(x), then F(k)=cG(k), cis a constant. 
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3.If f(x) = 





, then F(k)= 





Tels CONTIN 


A.If f(x) = g(h(x), then F(k)= S.GOH(k-D. 


1=0 
5. If f(x) =x", then F(k)=d(k—n), 6 is the Kronecker delta. 


n m m Àl 
md ras then F(k)=a> Ci(k+n r)! 


7.1f f(x) =a(1—bx) = 








(—b)’G(k +n—r), where aand 


bare constants 


3. Formulation of the Problem 


Free Vibration ofa Non-uniform Beam 


The governing differential equation for a non-uniform Euler beam shown in fig. 1 is 
given by: 








pA yey 4 Ža G: 2)- 0 (4) 


Where pis the density of the beam material, A(x)is the cross sectional area of the beam, 
v(x,t) is displacement of the beam, /(x) is the inertia at distance x from the left end of the 


beam and Fis young’s modulus of the beam. 





c ———_> 


Figure 1 Beam with variable cross-section 
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Boundary conditions: 


Case a. Simple -Simple Beam 








ô?v(0,t) d°v(L,t) 
0,t) =0; =0; L,t)=0; =0 
v(0,t) p v(L,t) a 
Case b. Clamped- Clamped Beam 
Op =0; OD o; gays; ZED o 
ox Ox 
Case c. Clamped- Roller Beam 
2 
v0,)=0; O20; wyso; ZEP o 
ð Ox 
Where Lis the beam length and fis the time 
4. Solution of the problem 
Assume that, the displacement of the beam is given by: 
v(x,t) = v(x) exp(iat) 
Where wis the natural frequency of the beam. 
Equation (4) can be conveniently written as: 
dV) yp dE) dW) pIE) dV) 
-pA(T)O? v(x) + E —— > 5 5 
dx* dx dx dx dx 











(5) 


(6) 


(7) 


(8) 


a0" (9) 
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With boundary conditions as, 


Case a. Simple -Simple Beam 








d’v(0) 














0)=0; =0; L)=0; =0 10 
v(0) : v(L) T (10) 
Case b. Clamped- Clamped Beam 
dv(0) dv(L) 
0)=0; =0; L)=0; =0 11 
v(0) = v(L) = (11) 
Case c. Clamped- Roller Beam 
dv(0) d°v(L) 
0) =0; =0; vw(L)=0; =0 12 
v(0) T v(L) TP (12) 


5. Dimensionless form 


Equations (9) can be conveniently written in terms of dimensionless variables as : 








d* y(x) 5 IS) yA) , d°S(x) PYA) m 
SO t e e ae OO (13) 
where 
EI (x) 3 V x 
= §(x)=(1- : a > 
El, (x) =(1— Bx) a e 
A(x) 


=S, SPs- px), 
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[í 
Q? = op Pe (the non-dimensional frequency of the beam), Ao, Ap, To are the 
0 


beam height , the cross sectional area and the inertia at the point at the left edge of the 


beam, and p is a constant equals zero for uniform beam. 


the non-dimensional boundary conditions are, 


Case a. Simple -Simple Beam 


d’y(0 d’y(1 
y(0) =0; 20 =0; Joso 22025 (14) 








Case b. Clamped- Clamped Beam 





y(0) =0; BOA: y(1) =0; A _9 (15) 
dx dx 
Case c. Clamped- Roller Beam 
dy(0 d’y(l 
geo, ag. Gaze. “OOo (16) 
dx dx 


6. Applying the Differential Transform Method 


In solving the problem, governing differential equations are converted to algebraic 
equations using DTM method which must be solved together with applied boundary 
conditions. Applying the differential transformation method (theorem (1-7)) to the non- 
dimensional governing equation (13) yield, 
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Yea 


j (k +1)(k +2)(k +3)(k +4) 
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| QY, — p’ [(k-DRMK +1)(k +2) + 6K (kK +DK + 2) + 6(k + DK +2)]¥,,, 
+2B[k(k +1)(k + 2k +3) + 3(k + Dk + 20K +3)] Y, 


(17) 


Then applying the Differential Transform Method to the non-dimensional boundary 


conditions equations (14)-(16) yield, 


Case a. Simple -Simple Beam 


The DT of Eqs. (14) is written as 


leads to Y[0]=0 


y (0)= Şa +1)(k+2)¥[k +2]x* = (D(2)¥[2]x° + (DOBI + 3)(4)¥[2] x7 


+(AOYBI + SOYA + (6)(7)Y[S}° 


leadsto Y[2]=0 


leads to > YIk] =0 


k=0 


Taking Y[l]=c 


y (= 5 (k +1)(k + 2)¥[k + 2]x* = (1)(2)Y[2]x° + (2)(3)Y[3]x' + B)(4)Y[2]x° 


k=0 


+(AOYBI + SOYA + (6)(7)Y[S}° 


leads to Skik —1)Y[k]=0 


k=0 


Taking Y[3]=d 


(18.a) 
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Equations (18.c) and (18.e) may be written as 


aa, bb, |c] |O T 
cc, dd,||d,| |0 C9) 


since cz and d7are not zero, for a non-trivial solution to exist the determinant of the matrix 
must be zero, i.e. 


a, x dd, —cc, x bb, =0 (20) 
where aaz, bbzare the coefficients of cz and dz in the equation (18.c) and cc;, dd; are the 


coefficients of cz and dzin the equation (18.e). The root of Eq. (20) is the solution for case a 
of the problem. 


Case b. Clamped- Clamped Beam 


The DT of Eqs. (15) is written as 
y(0)= i Y[k]x* =Y[O]x° + ¥ [lx + ¥[2]x* + ¥[3)x? + ¥[4]x* + Y[S]x° +.....=0 


leads to Y[0]= (21.a) 


y (0)= 3 (k + DY[k +1]x* = (DYMX? + (2)Y[2]x' + OYB]? + (HY[4]x 


k=0 
+(5)Y[5]x* + (6)Y[6]x° +.....=0 
leads to Y[l]=0 (21.b) 
y(1)= ri Y[k Olx? +Y[l]x' + Y[2]x? + ¥[3]x° + V[4]x* +Y[S]x° +.....=0 
leads to SYI] =0 (21.c) 


k=0 


Taking Y[2]=c (21.d) 
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y()= S (k +DY[k +1]x* = (DY[U]x’ +(DY[2]x" + (3Y [3] + (HY[4]x* 


k=0 


+(5)Y[5]x* + (6)Y[6]x° +... 0 


leads to X kY[k]=0 (21.e) 


k=0 


Taking Y[3]= d (21.f) 


Equations (21.c) and (21.e) may be written as 


aa, bb, ||c,| |O J 
cc, dd,||d,| |0 (22) 


since czand d?zare not zero, for a non-trivial solution to exist the determinant of the matrix 
must be zero, i.e. 


aa, x dd, —cc, x bb, =0 (23) 


where aaz, bbzare the coefficients of cz and dz in the equation (21.c) and ccz, ddz are the 
coefficients of czand dzin the equation (21.e). The root of Eq. (23) is the solution for case b 
of the problem. 


Case c. Clamped- Roller Beam 


The DT of Eqs. (16) is written as 


leads to Y[0]=0 (24.a) 


y (0) = y (kK+DY[k +1]x* =(DY[1]x° + (2)Y[2]x' + (3) Y[3]x? + (4 Y[4 1° 


k=0 
+(5)Y[5]x* + (6)Y[6]x° +.....=0 
leads to Y[l]=0 (24.b) 
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leads to Y YIK =0 (24.c) 
Taking ae c (24.d) 
y @)= y (k +1)(k + 2)Y[k + 2]x* = OAXDY[2]x? + (2)(3)Y[3]x! + GXAY[2]x 

S +(4)(S)¥[3]x° + (5)(6)¥[4]x* + (6OY +.....=0 
leads to 5 k(k -1)Y[k] =0 (24.e) 
Taking me j d (24.f) 


Equations (24.c) and (24.e) may be written as 


aa, bb, ||c,| |O 3s 
cc, dd,||d,| |0 >) 


since czand dare not zero, for a non-trivial solution to exist the determinant of the matrix 
must be zero, i.e. 

aa, x dd, — cc, x bb, =0 (26) 
where aas, bbgare the coefficients of cz and d3in the equation (24.c) and cc3, dd; are the 


coefficients of czand dzin the equation (24.e). The root of Eq. (26) is the solution for case c 
of the problem. 


7. Results and Discussion 


For a beam with simply supported ends Q? = 97.4091 leads to the first non-dimensional 
natural frequency as Q = 9.8696, which agrees with closed form value [5]. 


For a beam with clamped supported ends, Q? =500.564 leads to the first non-dimensional 
natural frequency as Q = 22.3733, which agrees with closed form value [5] 
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For a beam with one end clamped and the other end free, Q? = 237.721 leads to the first 
non-dimensional natural frequency as Q =15.4182, which agrees with closed form value 


[5]. 


Table 1 compares the accuracy of the first three non-dimensional frequencies of Simply 
supported uniform beams for different number of terms N. It is observed that the 
convergence speed increases with decreasing frequency order. i.e, the first frequency 1 
needs 25 terms to reach exact solution , while , third frequency 03 needs 45 terms. Table 2 
gives results of the first three non-dimensional frequencies of Euler beams with uniform 
cross section and different cases of boundary conditions. In table 2 exact values are also 
listed for direct comparison. It can be demonstrated that the differential transformation 
method is an efficient method in solving the vibrations of beams with good accuracy. Table 
3 gives results of the first three non-dimensional frequencies of simply supported non- 
uniform Euler beams with different values of beta £. It can be demonstrated that the first 
three frequencies of the simply non-uniform beam decreases with increasing 2 due to 
decreasing the beam cross section. Tables 4 and 5 gives results of the first three non- 
dimensional frequencies of clamped-clamped and clamped-roller non-uniform Euler 
beams with different values of beta Æ which give different cross sections. It is also 
observed that the first three frequencies of the clamped-clamped and clamped-roller non- 
uniform beam decreases with decreasing of the cross section. 


Table 1. The first three Non-dimensional frequencies of Simple-Simple uniform Euler 
Beams for different number of terms N 















































Qi Q2 Q3 
N=10 9.8902098156 - - 
N=15 9.8696683979 37.2824413197 - 
N=16 9.8696683979 37.2824413197 - 
N=17 9.8696020437 40.064696792 58.0531514621 
N=19 9.8696044699 39.4169139196 - 
N=20 9.8696044699 39.4169139196 - 
N=25 9.8696044011 39.4784501712 87.8912222720 
N=30 9.8696044011 39.4784176959 88.8107229014 
N=35 9.8696044011 39.4784176044 88.8264496362 
N=40 9.8696044011 39.4784176044 88.8264396509 
N=45 9.8696044011 39.4784176044 88.8264396098 
N=46 9.8696044011 39.4784176044 88.8264396098 
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Table 2. Comparison of results for free vibration of uniform Euler beams. 


The first three Non-dimensional frequencies of uniform Euler Beams: 












































. Q 1 Q 2 Q 3 
Uniform Beam 
Relative Relative Relative 
(6=0) DTM Exact] DTM Exact DTM Exact 
error % error % error % 
Simple-Simple 9.8696 9.8696 0.00 39.4784 | 39.4784 | 0.00 88.8264 88.8264 0.00 
Clamped- 
22.3733 | 22.3733 | 0.00 61.6729 | 61.6728 | 0.00 120.9031 | 120.9030 | 0.00 
Clamped 
Clamped-Roller 15.4182 | 15.4182 | 0.00 49.9649 | 49.9649 | 0.00 103.9982 | 103.9982 | 0.00 





-The first three Non-dimensional frequencies of Non-uniform Euler Beams for different 


values of £: 


Table 3. Simply supported Non-uniform Euler Beams 

















































































































El (x Q Q Q 
apy) : : 
EI,(x) B=0 | B=0.25 | B=0.5 | B=0 | B=0.25 | B=0.5 B=0 | B=0.25 | B=0.5 
Simple-Simple 9.8696 | 8.5772 | 7.1215 | 39.4784 | 34.4062 | 28.9519 | 88.8107 | 77.3785 | 64.9802 
Table 4. Clamped-Clamped Non- uniform Euler Beams 
EI (x Q Q Q 
(x) =(1- Bx) 1 2 3 
EI (x) B=0 | B=0.25 | B=0.5 B=0 | B=0.25 | B=0.5 B=0 B=0.25 | B=0.5 
Clamped-Clamped 22.3733 | 19.4836 | 16.3356 | 61.6729 | 53.6971 | 44.9817 | 120.9031 | 105.2123 | 88.1593 
Table 5. Clamped - Roller Non -uniform Euler Beams 
El (x Q Q Q 
<a py) : 
ElL (x) B=0 | B=0.25 | B=0.5 B=0 | B=0.25 | B=0.5 B=0 B=0.25 | B=0.5 
Clamped-Roller 15.4182 | 13.9524 | 12.3001 | 49.9649 | 44.0199 | 37.5276 | 103.9982 | 91.2744 | 77.1247 



































Figures (2-4) show the first three mode shapes of uniform Euler beams with 
different boundary conditions. From figures (2-4) It can be seen that the first three mode 
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shapes drawn using DTM is very agreement with the mode shaped drawn using exact 
methods. Figures (5-7) show, respectively, the first, the second and the third mode shapes 
of simply supported non-uniform Euler beams with different cross sections. Also figures 
(8-13) show, respectively, the first, the second and the third mode shapes of clamped- 
clamped and clamped-roller non-uniform Euler beams with different cross sections. From 
figures (5-13)). It can be seen that the value of beta (G=0, 0.25 and 0.5) has a significant 
effect on the mode shapes and the deflection values. Deflection of non-uniform beams 
increases as cross section decreases (f values increases). For Euler beam with uniform 


cross section (@=0). 














Deflection 
Deflection 




















Fig. 2 The first three mode shapes of Simple -Simple Fig. 3 The first three mode shapes of Clamped- 


uniform-beam Clamped uniform-beam 











Deflection 





Deflection 





ii 
Ka 


| = | 0. Beta=0 
erry Beta=0.25 
= 
Exact f — 


-0.07 ; ; f -0.3 
0 0.2 0.4 0.6 0.8 1 “o of 02 03 04 05 06 07 08 09 


X X 










































Fig. 4 The first three mode shapes of Clamped- Fig. 5 The first mode shapes for three different values 
Roller uniform beam of Beta of Simple -Simple non-uniform beam 
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Fig. 6 The Second mode shapes for three different Fig. 7 The Third mode shapes for three different 
values of Beta of Simple -Simple non-uniform beam values of Beta of Simple -Simple non-uniform beam 
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Fig. 8 The first mode shapes for three different Fig. 9 The Second mode shapes for three different 
values of Beta of Clamped-Clamped non-uniform values of Beta of Clamped-Clamped non-uniform 


beam beam 
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Fig. 10 The Third mode shapes for three different Fig. 11 The first mode shapes for three different 
values of Beta of Clamped-Clamped non-uniform values of Beta of Clamped- Roller non-uniform beam 


beam 
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Fig. 12 The Second mode shapes for three different Fig, 13 The Third mode shapes for three different 
values of Beta of Clamped- Roller non-uniform beam values of Beta of Clamped- Roller non-uniform beam 


8. Conclusion 


Based on the results presented, it can be demonstrated that the differential 
transformation method is an efficient method in solving the vibrations of beams with good 
accuracy using a few terms. The frequency of non-uniform beam decreases with decreasing 
of the cross section area and its inertia. Also it is observed that the deflection of the non- 
uniform Euler beam is increased as the cross section is decreased. 
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